The Mesoscopic Kondo Box: A Mean-Field Approach 
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We study the mesoscopic Kondo box, consisting of a quantum spin 1/2 interacting with a chaotic electronic 
bath as can be realized by a magnetic impurity coupled to electrons on a quantum dot, using a mean-field 
approach for the Kondo interaction. Its numerical efficiency allows us to analyze the Kondo temperature, the 
local magnetic susceptibility, and the conductance statistics for a large number of samples with energy levels 
obtained by random matrix theory. We see pronounced parity effects in the average values and in the probability 
distributions, depending on an even and odd electronic occupation of the quantum dot, respectively. These parity 
effects are directly accessible in experiments. 

PACS numbers: 73.23.-b, 71.27.+a, 72.15.Qm, 75.20.Hr 



I. INTRODUCTION 

Many-body effects have been a key interest in condensed 
matter physics for many decades. A prime example is the 
Kondo effect^ that, in its original context, refers to an in- 
crease of the resistance with decreasing temperature below 
the characteristic Kondo temperature Tk in metals contain- 
ing magnetic impurities. The significant progress in the fab- 
rication of mesoscopic and nanoscopic systems has lead to 
many alternative realizations of Kondo systems. A nice 
example is the so-called Kondo box,^ where a finite num- 
ber of electrons, confined in a quantum dot, is coupled to 
a single magnetic impurity. Their discrete energy spectrum 
introduces the mean level spacing A as new energy scale. 
It has been shown that physical quantities strongly deviate 
from the metallic behavior for A > Tk, i.e., when the size 
of the Kondo cloud screening the impurity would become 
larger than the size of the system. For example, they show 
parity effects, i.e., characteristic differences for an even and 
odd number of electrons in the system. Kondo physics 
in the presence of a chaotic dot geometry and disorder, re- 
spectively, has been studied within the Kondo disorder model 
(KDM)^^~^^ using Anderson's poor man's scaling approacbii 
to calculate the Kondo temperature. However, poor man's 
scaling cannot describe the strong coupling phase below Tk, 
and non-perturbative methods are needed, for example quan- 
tum Monte Carlo simulations (QMC),^^ numerical renormal- 
ization group ^^'^^ or the non-crossing approximation.^^ These 
methods are either numerically expensive or not reliable for 
low temperatures T <C Tk?'^ The QMC has been used by 
Kaul et alP to calculate the local magnetic susceptibility 
Xo for three different realizations of a chaotic system show- 
ing that mesoscopic fluctuations lead to significant deviations 
from bulk universality. However, parity effects and system- 
atic studies of mesoscopic fluctuations were neither consid- 
ered here nor in the KDM studies. 

In the present Letter we use a mean-field approach which 
had been developed initially for the single impurity Kondo 
model in a macroscopic metal,^"^ and later adapted to some 
particular mesoscopic Kondo systems. ^'^^ We first introduce 
the method in the framework of the chaotic Kondo box model. 



We then discuss the resulting probability distributions of Tk, 
Xo, and the conductance G. Their knowledge and the char- 
acteristic differences between even and odd electronic fillings 
make the comparison with experiments a realistic endeavor 
for the near future. 



II. MODEL 

We investigate the Kondo box model :^ An electronic bath, 
e.g. quantum dot, with discrete energy levels £u coupled to a 
local spin S = 1/2. The system is described by the Kondo 
Hamiltonian,^ 



Aa^l^ + S • So . 



(1) 



Operator (dia) creates (annihilates) a dot electron with 
level index / = 1, . . . , and spin component a =t, |. The 
chemical potential /i, related to an external gate voltage, fixes 
the number of dot electrons Nc, which effectively accounts 
for Coulomb blockade. In addition, the case of a fixed /i = 
rather than a fixed Nc will be considered. The local spin den- 
sity So of the dot electrons at the impurity position is given 
by 



4;^04,) /2, 



(2) 
(3) 
(4) 



where the local electron operator doa- is related to the dot elec- 
tron operators by 



doa = ^Aid, 



la ■ 



(5) 



The complex coefficients Ai correspond to the non-interacting 
wave function amplitude of level si at the impurity site 



(El \Ai 



1) and are assumed to be spin-independent. 



The quantum dot is thus characterized by N energy levels 
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£i, distributed from —W/2 to W/2, and the corresponding 
amplitudes Ai. In this Letter, our focus will be on chaotic 
quantum dots where the si are realized within random ma- 
trix theory^^ to be a Gaussian orthogonal ensemble unfolded 
to a constant density of states, and the intensities \Aif, ob- 
tained within the random wave modelr^ are Porter-Thomas 
distributed. For reasons of comparison we also introduce a 
system with equidistant ei and constant Ai that we will refer 
to as the clean system, cf., e.g., Fig.[T] 



III. MEAN-FIELD APPROXIMATION 
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For each given configuration, we treat the Kondo interac- 
tion in a mean-field approximation.^"^ The magnetic impurity 
spin is thus described in terms of fermionic operators /^^^ 
obeying the constraint + = 1- The Kondo inter- 
action in ([T]) then reads 



JkS • So = ^ ^ fafcT'dl^'doa -^Yl 4a^0a • (6) 



Jk 



The first term describes the spin flip processes, while the sec- 
ond term corresponds to a local potential scattering and is not 
further considered in the following. The mean-field treatment 
of the Kondo box Hamiltonian invokes two approximations, 
(i) replacing the quartic spin flip term by a quadratic term, 
r^^if^doa- + h.c), where the effective hybridization r is 
determined self-consistently by minimizing the free energy, 
and (ii) introducing a static Lagrange multiplier A, in order to 
fulfill (on average) the constraint of a single impurity spin. So 
the Hamiltonian ([T]) of the Kondo box in mean-field approxi- 
mation reads H = + H^, with 

= E [(^^ - ^) ^L^^- + riA^fldi^ + h.c.)] - XflU . 
I 

(7) 

Hereafter, since the mean-field approximation decouples the 
spin components, we suppress the spin index. The mean-field 
parameter r, the Lagrange multiplier A, and the chemical po- 
tential /i, satisfy the self-consistency equations 



r 

1 

2 

Nc 



(/tdo)=T^Ar^Gif(zc^n) , 

I iujn 

{Pf) = Tj2Gs{ioJn) , 



(8) 
(9) 
(10) 



I ioJn 



with the temperature T and the fermionic Matsubara frequen- 
cies ujn = (2n + 1) ttT. Here, the thermal average (• • • ) is 
expressed in terms of one-particle Green's functions, which 
are calculated from the mean-field Hamiltonian (|7]) using 



FIG. 1: (Color online) Kondo temperature Tr, for the clean system 
(Nc even o, Nc odd •) and average Kondo temperature T^^^^, for 
the chaotic system (Nc even 0, Nc odd /i = x), as functions 
of A/T^"^^. 



equations of motion. We find 
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Gff(zCJn) , 



(11) 

(12) 



\r?\Ai\' 

iuJn^ fi- Si 



-Gff {iuJn) 



(13) 



IV. RESULTS 



We now present the results starting with a discussion of T^, 
followed by an analysis of xo and G. In contrast to previ- 
ous studies of chaotic Kondo systems^^ we are able to distin- 
guish between an even and odd number Nq of dot electrons, 
cf. Eq. (fTOl) . We find pronounced parity effects in the prob- 
ability distributions of Tk, Xo and G. Solving only Eqs. ([5]) 
and (O, we can furthermore address the case of a non-fixed 
Nc, which corresponds to fixing /i = as done by Kaul et 
al. P All the numerical results presented here are obtained for 
N = 100 levels, and an electronic filling Nc = 101 and 102 
in the odd and even cases, respectively. We checked that N is 
large enough for band edge effects to not play a role. 



A. Kondo temperature 

Within the mean- field approximation, the Kondo temper- 
ature Tk characterizes a transition between a high tempera- 
ture phase, where the effective hybridization r, vanishes, and 
a low temperature phase with r ^ 0. More accurate meth- 
ods would rather describe a crossover, around Tk, between 
weakly and strongly coupled regimes. Despite the oversim- 
plified description of the T > Tk weakly coupled regime. 
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the mean-field approximation provides a good estimation for 
Tk, and a good description of the physical properties in the 
T < Tk strongly coupled regime.^ Using Eqs. ([8l[TTl[T2l), we 
find the well known Nagaoka-Suhl equation^ for Tk 

^ = E^tanhr£i^V (14) 
Jk ^ si-^ \2Tk J 

In the bulk limit A/W ^ 0, one obtains the relation^ 

T^^^^ ^ l.Uy^Wyl^ exp {-W/Jk) . (15) 

In the following, we use T^^^^, rather than Jk, as the refer- 
ence energy scale characterizing the Kondo coupling. 

The Kondo temperature probability distribution, p{Tk), is 
studied by solving Eq. ([T4b for various realizations. The result 
is mostly a regular, smooth function, which is depicted in the 
main part of Fig. [21 In addition to the data shown, there is 
a non- vanishing contribution p{Tk = 0) for Nq even that is 
analyzed in the inset of Fig. [21 Note that T^®^" is the average 
Kondo temperature computed from the regular part of p{Tk) 
only. 

Figure [H shows T^®^" for a chaotic system as a function of 
A/T^^^^, compared to the clean system, each for an even and 
odd electronic filHng of the dot. In the bulk limit (A /T^^^^ 
0), the Kondo temperatures for even and odd Nq coincide as 
expected, while there is a pronounced parity effect for increas- 
ing A/T^^^, i.e., smaller systems and/or decreasing Kondo 
interaction strength. For the clean system, Tk vanishes in the 
even case at the critical value A/T^^^^ ^ 6.28.^^ In all other 
cases and for larger values of A/T^^^^, the (average) Kondo 
temperature scales with the level spacing A, which becomes, 
in this limit, the only energy scale of the problem. 

The inset of Fig. [21 shows the fraction Pq of unscreened im- 
purities at T = as a function of A/T^^^^. We find that it 
is always possible to form the Kondo singlet in the odd case 
(no unscreened impurities), since the chemical potential coin- 
cides with an energy level. In the even case, the Kondo effect 
disappears, on average, with increasing mean level spacing A, 

The parity effects visible in T^®^" are even more dramatic 
in the probability distribution p (Tk) shown in Fig. [21 for 
^/2^buik ^ 2 and 100. Whereas for A/T^^^^ = 2 the av- 
erage Kondo temperatures are almost the same, there is al- 
ready a significant difference in p {Tk)- Most remarkable is 
a power law scaling, p (Tk) cx T^^, over several orders of 
magnitude in the odd case. Its existence was pointed out pre- 
viously ^"^"^^ for a related system, the Kondo disorder model 
(KDM). The KDM describes a Kondo impurity in a bulk elec- 
tronic bath for which disorder can be included continuously. 
The chaotic Kondo box that we consider here can in many re- 
spects be considered as a KDM with a non-tunable, i.e., fixed, 
disorder strength. Instead, another parameter can be tuned, 
A/T^^^, allowing for a continuous connection between the 
bulk, A/W = 0, and the mesoscopic, A ~ j^^uik j-ggi^ies 
Using Eq. ([T4l) we can show analytically that the power law 
scaling of p (Tk), with the constant exponent u = 1/2, is a 
direct consequence of the Porter-Thomas distribution of the 
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FIG. 2: (Color online) Kondo temperature distribution for a chaotic 
system with A/T^^^^ = 2 (soUd lines), 100 (dashed lines) for fixed 
Nc and /i = 0. For Nc odd we see a power law p{Tk) oc T^^, 
indicated by the straight line. Inset: probability for vanishing Tk in 
a chaotic system as function of A/T^""^^ for Nc even (0), Nc odd 
(♦), and/i = 0(x). 

intensities \Ai\'^. It persists for A > T^^^^ and is in agree- 
ment with Cornaglia et al.}^ In the case of even occupation, 
a linear probability distribution, p {Tk) ~ Tk was found for 
small Tk,^^ which is also consistent with our results. For val- 
ues A < T^u^^ the level structure is no longer important, so all 
the probability distributions are similar to a Gaussian around 

^mean 
-'-K 

B. Magnetic susceptibility 

Within the mean- field approximation, the Kondo spin static 
susceptibility 

Xo(T)^i^'^^dT(S(T)-S(0)), (16) 

reads 

Xo (T) = T ^ Gff (icon) Gff i-icon) . (17) 

ioJn 

In the decoupled phase, T > Tk, we recover a Curie law 
Xo {T) = 1/4T, which characterizes a free spin 1/2. In order 
to describe the screening in the Kondo phase T < Tk, we 
analyse the local effective moment T xo {T), which is equal 
to 1/4 for a free spin, and vanishes for a fully screened state. 
Figure [3] (upper panels) shows the local effective moment for 
fixed Nc as well as for /i = at A/T^^^^ = 2. The clean 
system (circles) shows parity effects as expected.- In the odd 
case, the impurity spin is fully screened at T = 0, while in 
the even case it remains unscreened or partially screened for 
all temperatures. In the chaotic system the average T xo only 
slightly differs from the clean case, since the screening of the 
impurity spin is closely related to the number Nc of dot elec- 
trons. More information provides the probability distribution 
p (T Xo), showing that in the odd case the impurity is screened 
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at T = for all configurations in the Kondo regime. For Nc 
even all values between and 1/4 are taken as T 0. This 
characteristic parity dependence in p {Txo) should be directly 
accessible in experiments. The /i = case contains both fea- 
tures: an increased probability for a screened impurity and a 
non- vanishing probability for having an unscreened impurity. 

C. Conductance 

Here, we consider the chaotic Kondo box connected to two 
leads, denoted by a = L^R. We study the tunneling con- 
ductance through the box in the linear response regime. The 
Hamiltonian of the full system is obtained from the Kondo 
box Hamiltonian ([T]) by 

H ^ H^^Skcl^^Ckaa^Yl (^r/c^a^/caa + h.C.) (18) 
kaa Ikaa 

where the second term describes the c— electrons of the leads 
and the third term the dot-lead coupling. The leads are as- 
sumed to be ideal metals, therefore, the tunneling couplings 
do not depend on the momentum k, tfj^ = tf. In order to 
mimic the chaotic nature of the Kondo box, we also assume 
that the tunneling couplings are randomly distributed, with 
second moments {tftf, ) = t'^Sw (here, S denotes the Kro- 
necker symbol, and t is a characteristic tunneling energy). In 
the presence of a finite bias voltage between the leads V, the 
current through the Kondo box is2^ 

J{V) = ^ j duj [fn {u) - fL {co)] Im [Tr {TG}] , (19) 

where /r/l = / (cj ± eF/2) is the Fermi function in the 
leads, and G is the Green's function matrix of the dot elec- 
trons, in the presence of the leads. 

r^, =27r^p(efc)«, (20) 

k 

is the tunneling matrix, where p denotes the density of states 
of the leads. Hereafter, we approximate the tunneling matrix 
by its average value, T^, ^ jS^, where 7 ~ pt^. A fully 
self-consistent mean-field treatment of the system with a finite 
bias voltage would require to take into account the renormal- 
ization of the mean-field parameters, in the presence of the 
leads, similarly to the approach of Aguado et ah P Here, we 
consider the linear response regime, V 0, as well as the 
tunneling limit, 7^0. The mean-field parameters r, A and 
/i are thus not modified by the leads. Furthermore, within the 
mean-field approximation, the effective Kondo box Hamilto- 
nian d?]) is non-interacting. Therefore, the Dyson equation for 



G reads 

G-^ = Gq^ -ir, (21) 

where the Kondo box (i.e. without the leads) (i— electron 
Green's function Go is given by Eq. ([T3l) . 

The conductance G = dJ/dV is shown in Fig. [3] (lower 
panels) as a function of T/Tk- There are clear parity effects: 
for Nc even, the low temperature limit is G = 1 (in units of 
the quantum of conductance 2e^//i), no matter what the un- 
derlying realization is. However, for an odd occupation num- 
ber the conductance goes to G = 0, since the Kondo singlet 
blocks the transport channel. The low temperature limit for 
/i = depends on the distance between the chemical poten- 
tial /i and the closest energy level. 



V. CONCLUSIONS 

We have presented a mean-field approach to the mesoscopic 
Kondo box problem that allows a very efficient computation 
of all physical quantities of interest - Tk, Xo, and G - that 
are easily, and on an equal footing, accessed via the one-body 
Green's functions. In contrast to other methods, we are able 
(i) to reach very low temperatures T ~ 0.001 Tk within rea- 
sonable computation time and (ii) to calculate the probability 
distributions based on a large number (at least 50 000) of re- 
alizations of the chaotic Kondo box. Our results agree with 
those of other approaches when available. Concretely, we 
find deviations from the bulk system to occur in form of pro- 
nounced parity effects. For realizations with an odd number 
of dot electrons, we confirm and refine the power law distribu- 
tion of p (Tk)- The significant parity effects in the magnetic 
susceptibility and the conductance provide the basis for a di- 
rect comparison with experiments. For example, the expected 
spread in the experimental values in quantum dots with a cer- 
tain (even or odd) number of electrons, fixed via Coulomb 
blockade, will be either large or small, and is therefore already 
accessible by measuring few samples. 
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